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Abstract: Zamolodchikov's famous analysis of the RG trajectory connecting succes- 
sive minimal CFT models Mp and Mp_i for p ^ 1, is improved by including second 
order in coupling constant corrections. This allows to compute IR quantities with next 
to leading order accuracy of the 1/p expansion. We compute in particular, the beta 
function and the anomalous dimensions for certain classes of fields. As a result we 
are able to identify with a greater accuracy the IR limit of these fields with certain 
linear combination of the IR theory Mp_i. We discuss the relation of these results with 
Gaotto's recent RG domain wall proposal. 
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Introduction 

In his famous paper [1] A. Zamolodchikov has investigated Renormalization Group 
(RG) flow from the minimal model Mp to the Mp_i for large p ^ 1 caused by the 
relevant operator 0i 3. Two main circumstances made it possible to investigate this 
RG flow using a single coupling constant perturbation theory. First, the conformal 

dimension of the field 0i^3, Ai 3 = 1 = 1 — e (see Appendix A) is nearly marginal 

when p ^ 1 and second, the Operator Product Expansion (OPE) of this field with 
itself produces no relevant field besides the initial one and the unit operator. The 
method of A. Zamolodchikov not only allowed to identify the IR theory with Mp_i, but 
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also provided detailed description how several classes of local fields behave along the 
RG trajectory. The analogous RG flow for the N = 1 super minimal models has been 
investigated in [2]. 

The main purpose of this paper is a sharpening of Zamolodchikov's analysis, by 
inclusion of the second orc^er perturbative corrections. It is interesting to note that in all 
cases we have investigated, the rotation matrix (in the space of flelds), that diagonalizes 
the matrix of anomalous dimensions, does not receive 1/p or l/p"^ corrections. So an 
interesting question arises, if any higher order corrections appear at all. 

As intermediate results, in this paper we have found several four-point correlation 
functions in large p limit (see formulae (C.l)). 

The initial motivation to carry out these computations came from the recent ap- 
proach to this RG flow by D. Gaiotto [3]. Using Goddard-Kent-Olive construction, 
Gaiotto has constructed a non-trivial conformal interface between two successive mini- 
mal models and made a striking conjecture, that this interface is the exact RG domain 
wall which encodes the map between the UV and IR flelds. A stringent test in favour 
of this conjecture provided the comparison in leading order of the 1/p expansion of the 
mixing amplitudes deduced from the RG domain wall approach with those obtained by 
Zamolodchikov using flrst order perturbation theory. In this paper we show that this 
mixing coefficients, computed with the help of the perturbation theory up to second or- 
der, unlike those obtained from the Gaiotto's conjecture, do not receive any corrections 
up to 1/p^ order. Nevertheless, this discrepancy might be attributed to the renormal- 
ization scheme which is adopted here following Zamolodchikov. Presently the author 
of this paper does not have any clue how to take into account possible dependencies 
on the renormalization schemes in order to be able to make any conclusive statement 
about Gaiotto's conjecture beyond the leading order. 

The paper is organized as follows. 

In Section 1, we develope some technical tools, necessary to carry out second order 
in coupling constant calculations. 

In Section 2 the /3-function and Zamolodchikov's c-function [4] are computed with 
next to leading order accuracy. The critical value of the renormalized coupling constant, 
the slope of the /3-function as well as the c-function at this point are calculated. The 
results of these computations confirm that the second order contributions too perfectly 
match with Zamolodchikoved's conclusion that the IR fixed point corresponds to the 
CFT Mp_i and that the UV field 0i^3 fiows to the field 03^i of the IR theory. 

Section 3 is devoted to the renormalization of several series of local fields and 
calculation of their anomalous dimensions. 

In Section 3.1 we investigate the renormalization of the fields 0n,n- 
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In Section 3.2 the renormalization of the fields 0n,n+i and 0n,n-i is discussed and the 
matrix of anomalous dimensions is found. At the fixed point the matrix of anomalous 
dimensions is diagonalized and its eigenvalues are calculated. 



In all cases the predictions of Zamolodchikov successfully withstand the next to 
leading order test. 

In Appendix A some basic facts about the minimal models of 2d CFT are reviewed. 
The Appendices B and C are devoted to computation of the integrals used in the main 
text. The Appendix D comments how to calculate the large p limit of those four point 
correlation functions used in the main text. 

1 Perturbation theory in second order 

Suppose the (Euchdean) action density is given by 



with Tio being the UV CFT action density, a relevant local spinless field and A the 
coupling constant. Then for a two-point function up to second order we'll have 



In this paper we consider a theory, whose UV limit is given by the minimal CFT model 
Mp with p ^ 1 and the perturbing field (p = (pi^^. Leading order corrections in this 
theory has been investigated by A. Zamolodchikov [1]. Second order computations are 
more complicated. Indeed, not only the knowledge of four point correlation functions 
which, in general are quite non-trivial in a CFT [5], but also integration over two 
insertion points is required. Fortunately, as we demonstrate below, the conformal 
invariance allows to perform integration over one of the insertion points explicitly. 
First let us notice that translational and scale invariance can be exploited to locate the 
points yi, y2 at = 1 and ?/2 = without loss of generality: 



In Section 3.3 the same steps are performed for the fields 0n,n+2, dd<Pn,n and 0„,n+2- 



U{x) = Hq{x) + A0(a;) 



(1.1) 



(0l(yi)02(y2))A 




(1.2) 




(1.3) 
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(here and below I frequently use the shorthand notation Xi2 = Xi — X2, yu = Hi ~ 1/2 
et.al.). Any four-point function of primary fields in a CFT essentially depends only on 
the cross ratio x = f^^f^ of the insertion points and its conjugate [5] 

{MX1)MX2)MX3)MX4)) = {xuXur^^'{x24X24)'''^^'-^'-''' (1.4) 

X (x3#34)^^+^^-^^'^Ha;23X23)^^-^^-'^^-'^^G(x,x), 

where it is assumed that the fields are spin-less (i.e. = Aj). Specifying the insertion 
points as xi = x, X2 = 0, X3 = 1 and X4 = 00 we get 

Gix,x)= lim (X4X4)'^M01W02(O)03(1)04(X4)) = (01(X)02(O)03(1)04(OO)) (1.5) 

Alternatively specifying xi = 1/x, X2 = oo, 0:3 = 1 and X4 = and comparing with 
(1.5) we get the identity 

{Mx)Mo)M^)Moo)) = {xx)-^^'{M^/x)Mo)M^)Moo)) (i.e) 

which is useful when investigating the correlation functions at large x. After application 
of (1.4), (1.5) to the four-point function ((/)(a;i)(/)2(O)0i(l)(/)(a;2))o and introduction of 
the new integration variables 

Xi(l-X2) . 
> Xi, 1 — X2 -> X2- 

X12 

the integrals on r.h.s. of the eq. (1.3) become partly disentangled 



(0i(l)02(O)0(xi)0(x2))orf'a;id^X2= / /(a;i)(0(xi)02(O)0i(l)0(oo))orf'a;i (1.7) 
where 

Hx) = J {yyr-\{l -y){i- y)f-\{x -y){x- y)r<fy, (1.8) 

with parameters 

a = ei2 + 2e; 6 = £21 + 2e; c = -2e (1.9) 
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where e = 1 — A, t\^2 = 1 ~ ^\,2 are the complementary dimensions. Fortunately the 
integral (1.8) can be expressed in terms of hyper-geometric functions (see appendix B): 

^ ![7(M^±4|F(l-a-6-c,-c,l-a-c,x)|2 (1.10) 
7(a + + c) 

^7(1 + c)7(a) 1 _ 6, 1 + a + c, a:)p 

7(1 + a + c) ' ^ ' ' ^' 

7r7(a)7(6 + c)|^,^ „ ^ „ „ ^ „ i ^m2 



7(0 + 6 + c) 



|F(1 — a — 6 — c, — c, 1 — 6 — c, 1 — a; ) 



7(1 + 6 + c) 

7r7(a)7(6),^^^,^ _ , ^ ^ ,^ 2 



7(a + &) 



|x'^F(a, — c, a + 6, 1/a;) 



+ -l(l±fM^±^|^»+^+e-i^(, - a - 6 - c, 1 - 6, 2 - a - 6, l/x)p 
7(a + 6 + c) ^ ' 

where 7(x) = r(x)/r(l — a;) and F{a,b,c,x) is the Gauss' hypergeometric function. 
Above three expressions for I{x) are convenient when exploring the regions x ~ 0, 
X ~ 1 and X ~ 00 respectively. Note also that these expressions make explicit the 
single- valuedness of /(x). Specifying the choice of parameters to (1.9) and applying 
the identity 

F{a, b, c, x) = (1 - zy-''-^F{c -a,c-b,c, x) 
to the second term of the second equality, the eqs. (1.10) can be rewritten as 

vr7(2e + e2i)7(ei2) I j;.., 001, ^2 

= — |F(1 -2e,2e, l + e2i,x)| 

7(2e) 

7(2e) 

7r7(2e ei2)7(e2i) I p,.-, 001, 1 m2 



7(2e) 



■|F(1 -2e,2e,l + ei2,l-x) 



+ ""^'^V;f^^"^ |(x/(1 - x))-F(2e, 1 - 2e, 1 + ^2. 1 - x)^ 
7(2e) 

- "^^'^ + t^|x--F(2e + e,2,2e,4e, l/x)^ 



7(4e) 

+ ^^^%7T^I^''''^(1 - 2e, 1 - 2e + ei2, 2 - 4e, l/x)p (1.11) 
7 (2e) 



It is worth noting that in the case when ei2 = ei — €2 = only the third expres- 
sion is manifestly nonsingular, the previous two expressions require a subtle limiting 
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procedure. Thus it is better to employ the third expression: 

7172 (2e) 



I{x) 



+ 



7(46) 



\x 
v. 



7^(2e) 



2^F(2e,2e,4e,l/x)|2 
\x'^''^F{l - 2e, 1 - 2e, 2 - 4e, 



:i.l2) 



Let us investigate the behaviour of (1.12) at x ~ 1. Using standard formulae for the 
analytic continuation of the hypergeometric function with parameters satisfying the 
condition a + b — c E Z (see e.g. [7]) one can get convinced that 



I{x) ^ 7r(a; + ^ - 4) + TT (1 + 2e(2e - l){x + x - 2)) 
X {2-\og\y- Ip - 27rcot(27re) - 4V^(2e) -47) 



:i.l3) 



where 7 = 0.577216 ■ ■ ■ is the Euler constant and the omitted terms are at most of 
order Ix — Ip log |x — 1| in a; — 1 limit. There is no need to investigate the limit a; — ?■ 
separately since the obvious symmetry of I{x) with respect to s ~ 

immediately ensures 



I{x) ^ -7r(2 + X + x) + TT (1 + 2e(l - 2e)(x + x)) 
X (2 - log |xp - 27r cot(27re) - 4^/'(2e) - 47) 



:i.l4) 



2 /^-function 



In this section we calculate /^-function up to 1/p^ ~ corrections for the small values 
of the (renormalized) coupling constant (of order e or smaller). As it will become 
quite clear later, for this purpose one should evaluate the integral (1.7) for the case 
01 = 02 = and I{y) given by (1.12), with the accuracy ~ 1/e. Our strategy will 
be as follows: separate in the integration region the discs Di^q = {x G C| \x\ < I}, 
A,i = {x G C I |x — 1| < /} and -Dz,oo = {x G C | |x| > 1//} where / is an intermediate 
length scale such that < /q ^ exp(— 1/e) ^ / ^ 1 /q being the ultraviolet scale. For 
the integral outside this discs we will safely use the small e limits of the correlation 
functions given in the appendix while inside the discs we'll explore (exact in e) OPE. 
We will see that the trace of the intermediate scale / will be washed out entirely from 
the final result. In present case the e = limit of the four-point function is given by 
(see appendix C) 



(0(x)0(O)0(l)0(oo)) 
1 - 2x + 3x2 - 2a;3 ^ 
x^fl — x)2 



(2.1) 



16 


1 3x < 2 

-L 2 ~'~ 4 


2 


5 


X 


y 


x(l — x)2 


H 


' 9 


1 — X 
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With required accuracy I{x) ~ 7r/e. It is convenient to carry out the integration 
in radial coordinates x = rexp{iip), x = rexp{—i(p), d?x = rdrdip. The resuh of 
integration over the angular variable ip will depend on the region, where the radial 
coordinates takes its value 



ReSx=oR{x, r^/ x)/x + ReSx=r^R{x, r^/ x)/x, if r < 1 
ReSx=oR{x, r"^ /x)/x + ReSx=iR{x, r^/x)/x, if r > 1 



for arbitrary rational function R{x,x) with poles located at x = or x = 1. In 
particular when R{x,x) is the r.h.s. of the eq. (2.1) we get 



3^10 _ 9^8 ^ 25r^ - 23r^ + 7r^ + 3 

3r4(l -r2)3 
3^10 ^ _ 23r^ + 25r^ - 9r^ + 3 



if r < 1 

if r>l (2.3) 



3^4(^2 _ lys 

After performing the remaining elementary integration over r we finally get 

Ott tt 337r 327r 
(0(x)0(O)0(l)0(oo))rf2x =- + — -- \og{2y') + ■■■ (2.4) 



where (see Fig. 1 ) 

^1,10 = (I^l-/„,o\A,o) U (Di//,o\i^l+Zo,o) 

and dots stand for the negligible terms of order / or Iq/I. There is a subtlety to be 
treated carefully here. The fact that the part of the white narrow ring (of width 21q) 
outside of the red circle is missing from the integration region is insignificant 
since its inclusion would produce only negligible terms of order Iq. Instead we have to 
subtract the contribution of two lens-like regions of fli^i^ included in the red circle (as 
already said, the contribution coming from the regions around the singular points will 
be computed separately exploring OPE), 
a) Integration over lens-like regions 
Expanding (2.1) around x ~ 1 we get 

1 2 2 16 , , 



|X-1|4 (X-1)2 (X-1)2 3|X 

where only the singular terms whose integrals over the regions around x = 1 diverge 
are presented. The integrals of such terms have been evaluated in Appendix D. As a 
result the contribution of the lens-like regions, to be subtracted from the r.h.s. of eq. 
(2.4), is equal to 
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Figure 1. i^n^ is the gray region 
b) Contributions of the regions around singularities 

It remains to calculate the the contributions of the regions Di^, Di^i, Di^^o to the 
integral (1.7). Evidently the first two give identical contributions, so let's consider the 
region Di q for definiteness. To calculate the four-point function in this region we apply 
the OPE (all the structure constants we use in this paper can be extracted from the 
general formula (A. 5)) 

0(x)0(O) = (xxr'^il + ■■■) + 4'3)(i,3)(^^)"'' (</'(0) + ■ ■ ■ ) (2-7) 
Taking into account (1.14) and that 



c£3)(i,3)C(i,3)(i,3)(i,3) = y (1 - 3e + 0{e')) {21 



one easily gets 



l{x){(j){x)(j){0)(f){l)(f){oo))d^x 

, 327r^(log(0-3) 327r^ 

where the first two terms come from the identity and last two terms from the (p field 
channel of OPE (2.7) respectively. 

c) Contribution o/x ~ cxd 
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At large x first we make use of eq. (1.6) to pass to tlie inverse variable 1/x ~ and 
only afterwards we apply OPE. The calculation is similar to the previous case, the main 
difference being the fact that at this limit I{x) becomes simply -(xx)~^^. The result is 



J(x)(0(x)0(O)0(l)0(oo))c/'x 
vr^ vr^ 16n^ { I 3 - 21og/ 



(2.10) 



Lt's pick up all the ingredients: (2.4) times I{x) ~ minus (2.6), plus twice (2.9), 
and plus (2.10). We get 

— + ^ + 0(e°) (2.11) 
As promised, the / dependence disappeared. Presence of the divergent term -|?^ also is 

'o 

not surprising. In our naive regularization scheme we could cancel this infinity adding 
an appropriate, proportional to the area counter-term in action. In fact, if we would 
have been able to treat the integral (1.7) analytically as continuation from a region 
of parameters where integral converges, then this kind of non-analytic divergent term 
couldn't emerge at all. In what follows we'll simply drop out such terms without further 
ado. 

Thus for the two point function G{x, A) = {(j){x)(j){Q))\ we get (c.f. (1.2)) 

G{x, A) = (xx)-2+2^ (^-\^ ~ ^ + ^^^)) 

A^ /SOvr^ SStt^ \ \ 

+ y — + ) ) (^^) +■■■) (2-12) 

Following A. Zamolodchikov let us introduce a new coordinate g ("renormalized" cou- 
pling constant) in the space of one-parameter family of theories (1.1) instead of the 
initial coupling A and introduce the local field 0^^^ = dgl-L (according to (1.1) the initial 
"bare" perturbing field (j) = dyhi). The new coupling g is fixed by the requirement that 
the two point function G{x,g) = {(f)^^\x)(f)^^\0))x satisfies the normalization condition 

G{l,g) = l (2.13) 

Then the /3-function can be computed from the identity (see [1]) 

e(x) = eX(j){x) = I3{g)(f)^'^\x) (2.14) 
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where G is the trace of the energy-momentum tensor. Combining (2.13) and (2.14) one 
easily finds 

d)^g=^G{l,\) (2.15) 

and 

/?((?) = eX^/W^) (2.16) 

The equation (2.15) allows one to express g in terms of A (the integration constant can 
be set to zero so that the unperturbed CFT will corresponds to (? = 0) 

^ = A - ^ f - - 3 + 0{e)\ + ^ f 4 - - + 0(1)1 + O(A^) (2.17) 



or, inversely 

A = ^ + ^ (' - 3 + 0(e)) + ^ (I - ' + 0(1)) + 0{g^) (2.18) 
Inverting and replacing in (2.16) A in favour of g we get 

m = - ^ (2 - 36 + 0(6^)) - ^ (1 + 0(6)) + ■ ■ ■ (2.19) 
The equation 

I5{g*) = (2.20) 

admits a nonzero solution 

27r/ = 736+ ^6^ + 0(63) (2.21) 

so we have a non-trivial infrared fixed point. In [1] this fixed point has been identified 
with the minimal model Mp_i and the local field (f)^^'^ with the field (f)fi Now we are 
in a position to check this identification more accurately. The anomalous dimension of 
0^^*) is related to the slope of /3-function 

A* = 1 - dg(3{g) = 1 + 6 + 6^ + 0(6^^) (2.22) 

which matches to the conformal dimension of computed from the Kac formula 

(A. 2). Also the shift of the central charge [4] 

c*-c, = -12n' (3(g)dg = -—- — + 0ie') (2.23) 
neatly matches to the exact expression 

12 _ 36^ 

" "(2-6)(l-6)- 
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3 Field renormalization and and the UV - IR map 



In this section we calculate the matrices of anomalous dimensions for several classes of 
fields. Diagonalization of these matrices at the IR fixed point provides a detailed map 
between the UV local fields and their image under RG flow in the IR theory. 

3.1 Primary fields 0„ „ 

This is the simplest case to analyze since the fields never get mixed with other 
fields [1]. This follows from the structure of the OPE involving the perturbing field 0i,3. 
The subspace of fields which is generated by the field and is closed w.r.t. OPE with 
01^3, doesn't contain any other field with a dimension close to that of A„ „ = O(e^). 
We are going to calculate corrections to the anomalous dimension up to the order e^. 
That is why for the present purpose the knowledge of the four point function 
(0(x)0n,n(O)0n,n(l)0(oo)) up to correctiou is required. As in previous case, to find 
this correlation function we first used AGT relation to find the relevant conformal 
blocks up to sufficiently large level (actually the the computations were performed to 
the order terms). Expanding a conformal block up to and examining first few 
coefficients of the resulting power series in x it is possible to guess the entire power series 
and identify it with some elementary function. Having in our disposal the expression 
for the correlation function we then checked that it satisfies all the nontrivial physical 
requirements: the single-valuedness and the compatibility with OPE around the points 
X ~ 1 and X ~ oo. Here is the final expression (see Appendix C) 



(0(x)0„,n(O)0„,„(l)0(oo)) (3.1) 

= 1+ .. +41og^ +0e3 

12 \2x[x — 1) 2x[x — 1) X J 

From eq. (1.12), up to order e, I{x) is equal to 

I{x) = ^-27rlog|(l-x)x|+87relog|a;|log|l-x|+0(e^) (3.2) 

Now we are ready to perform integration over the region Qi^^i (see Fig. 1). Since the 
singularities at x ~ and x ~ 1 are integrable, we can put Iq = 0. As in Section 2 
the integration over the angular variable should be performed separately for the cases 
< |x| < 1 and |x| > 1. Integration of rational expressions we have already discussed 
earlier. As about the logarithmic terms, they can be easily handled first expanding 
into power series in x when |x| < 1 or in 1/x otherwise. Then we pass to the radial 
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integration. Both steps are elementary and we present only the final result 

/(a;)(0(a;)0„,„(O)0„,„(l)0(oo))ci'x ^ ^— f logy + ij (3.3) 



'0 



+ - (2 + 4e + (4 + 8e) log/ + Selog^ /) - + 4e) ^— ^ + — 

Due to the already mentioned mildness of singularities at and 1 the only remaining 
contribution to be taken into account comes from the neighbourhood of oo i.e. from 

,oo • 

At large x it is convenient to employ eq. (1.6) 

(<^(x)0„,„(O)(/)„,„(l)0(oo)) = (xx)-2^(0(l/a;)(/)(O)0„,„(l)(^„,„(oo)) (3.4) 

and apply the OPE (2.7) with x replaced by 1/x. The correlation function decomposes 
into a sum of two partial amplitudes one corresponding to the identity and the other 
to the field (p. 
a) Contribution of identity 

The prefactor (xx)"^"^ in (3.4) compensates the factor (xx)^^ accompanying the iden- 
tity operator in OPE and with sufficient accuracy we can replace this partial amplitude 
by 1. It is straightforward to expand I{x) given by (1.12) at large x keeping only those 
terms which after integration may produce non-vanishing terms in small / limit 

^7'(2e). ._2. A , ^\ , ^\ , ^7(4e-l), .2e-i 



Integrating this expression over the region -D«,oo\-Dio,oo, dropping out, as earlier, all 
singular in /q terms and expanding the result up to the linear in e terms we get 

vr^ ^ 27r^(21og(/) + l) , 47r^6 (/^ - 2 log^(/) - 2 log(/) - l) 
-^ + ^ P + P (3.6) 



b) Contribution of the field 0i,3 
This contribution is 



e 



where ^{xx) is just the function I{x) with required accuracy, (xx) ^"""^^ is the pref- 
actor of (3.4), (xa;)^""^ comes from OPE and the squared structure constant is equal 
to 



'-^{l,3)(n,n)(n,n) " q <~ U{e ) {6. 
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The integral is converging at the hmit Iq ^ 0, so, we may perform integration over the 
entire region -D;,oo- The result reads 

7r2(n2-l)(l-e + 2elog(/)) 



(3.9) 



6 

The sum of all contributions (3.3), (3.6), and (3.9) is 

TT^ (^2 - 1) (2 + 5e) , , 

— ^ '- + 0{e') (3.10) 

Combining this with the first order in coupling constant contribution 

{<PnA^)<Pn,nm{^))d'x = ^ " 1) (2 + Se) 6 ^ ^^^^^^ 

8v 3 



C(l,3){n,n){n,n) = /- ^ (3-12) 



where the value 

(n2 - 1) (2 + 5e) 
for the structure constant is inserted, we get 

Gn{x,X) = (0n,n(a;)0„,n(O))A 

/ f, X ^(^'-l)(2 + 5e + 0(e2))e 

= (xx) ^^"'^ 1 - A ^ — ^ (xxY 

V 8^3 

^^ .^(„^-l)(2 + 5e + 0(a) (^^)..^...^ (3,13) 

Let's introduce the renormalized field (pnji = B{X)(f)n,n requiring that the two point 
function Gn{x,g) = {4>n}n{x)(l)n}i{0)) x satisfies the normalization condition 

Gn(l,^7) = l (3.14) 

so that 

Then for the anomalous dimension we get 

Alfl = A„,„ + eA dx log B = A„,„ - ^ 9aG„(1, A) (3.16) 
In view of (2.18) we find 

A£ = A„,„ (3.17) 
7r^(n2- 1)6^2 + 56 + 0(a) _ jn' - 1) e\l + O je)) , 
16^3 8 ^ 
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So that at the fixed point 



64 ^'-''^ 

which completely agrees with the dimension An,n^^ of the field (p^^n^'^ in the minimal 
model Mp_i. Thus the conclusion of A. Zamolodchikov that under RG the UV field 
(pnjt fiows to IR 4'n,n^^ is robust also against our second order test. 

3.2 Renormalization of the fields 0„,„+i and 0„,n-i 

Already in this case we encounter the phenomenon of mixing. The OPE 0i,30n.n+i 
produces besides 0„,n+i also the primary field 0n,n-i both having dimensions close to 
1/4 in large p limit. Thus we have to consider the correlation functions 
(0(x)0„^„-i-i(O)0„^„-i-i(l)0(oo)) with all four possible choices of signs. The strategy is 
exactly the same as in previous sections and for each choice we will follow the steps 
performed in Section 3.1. 



3.2.1 Correlation function (0„^„_|_i(1)0„^„+i(O))a 

a) Contribution of the region Qi^^^ 
This is given by the integral 

I{x){(f){x)(f)n,n+liO)(pn,n+li^)(pioo))£x 



(3.19) 



The large p limit of the four-point function found from AGT relation is (See Appendix 
C): 



(0(x)0„,n+i(O)0„,„+i(l)0(oo)) 



2{n + 2) 
3n 



X 



xix — V, 



9 1 

>Xj I n 



x{x — 1) 



+ 0(f 



(3.20) 



With required accuracy l{x) can be replaced by ^. The integral (3.19) can be performed 
using the technique already explored in computing (2.4) or (3.3). The result is 



.10 



TT 



/(x) (0(x)0„,„+i (O)0n,n+1 (l)0(oo))(i^X 

7r(13n + 20) log(/) 7r(5n + 4) log (2/o) 



6n 



(3.21) 



6n 

b) Contribution of lens-like regions 

Near x ~ 1 the r.h.c. of eq. (3.20) becomes i2n\x'-i\'2 ' hence the contribution of the 
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lens-like regions near x ~ 1 to be subtracted from the r.li.s. of eq. (3.21) is equal to 
(see (D.9)) 

''^""^'\oJ±) (3,22) 



6ne \2^o, 
c) Contributions of the regions Di q, Di i and Di ^o 

Contributions of Di q and Di^i obviously are identical so we will concentrate on Di^o 
only. The relevant OPE is 

+ ^a57ni+l)(^^)^"'"-^-^"'""^"^ (0n,n-l(O) + ■ ■ ■ ) (3.23) 

Taking into account that in this region I{x) ~ 7 — 7rlog(|xp) (see (1.14)) and that 
r2 _ (^ + 2y (1 - (2n - 1)6) 

<-^(l,3)(n,n+l)(n,n+l) " ^ ) 

C - + , oie^\ r3 241 

'-^{l,3)(n,n+l)(n,n-l) " '^^^ T U{e ) l<J-^^j 

we get 

/ X , / X / X / X / X V ,9 TT^fra + 2)^ / 1 1 — 2r;, + log / 



L 



, 27r^(n^-l) 1^ 72 + 4 ^ n + 4 + (n + 2)Mog^ ^^ 
+ 3n2(n + 2)2 V e 7 ^ ^ 

Above two terms come from two primaries 0n±i appearing on the r.h.s. of the OPE 
(3.23). 

The contribution from the region Di ^o is completely analogous to the case of cor- 
relation function (00) a discussed in Section 2. The only difference is that now the 
contribution of the field which appears in u-channel OPE is proportional to 

^(1,3) ^ _2(n + 2HW!i±lM , nr.2N ^.a^ 

'-^(l,3){l,3)'-^(1.3){n,n+l){n,n+l) " 1^ H U\e ) 

instead of (^(^i 3)(i 3)(i 3) ~ ^^(^-3*^) _ xhe result is (c.f. eq. (2.10)) 

/(a;)(0(a;)0„,„+i(O)0„,„+i(l)0(oo))rf^x 
^ vr^ , 27r^(n + 2) (\ n + l-21og/\ 

^ eP^ 3n [e^ e J ^^"^^^ 

Picking up all the contributions: (3.21) minus (3.22) plus twice (3.25) and plus (3.27), 

we get 

TT^ (3n^ + 24n2 + 64n + 44) An'^in + 1) {n^ + 7n^ + Un + 5) ^ 
Snln + Ij^e"^ Sn[n + 2j^e 
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3.2.2 Correlation function (0n,n-i(l)0n,n+i(O))A 

a) Contribution of the region Qi^i^^ 
is given by tlie integral 



/(x)(0(x)0„,„+i(O)0,,„_i(l)0(cx)))d'a; (3.29) 
The large p limit of the four-point function now is (see Appendix C): 



(0(x)0 (O)0„,„_i(l)0(cX))) 



Vn"^ - 1 2x-l ^ 



3n 



x{x — 1) 



+ 0(e) (3.30) 



Using the last equality in (1.11) we see that I{x) ~ (^iq^^'^^^ result of the integral 

(3.29) is 

l{x){(t){x)(f)n,n+l{0)(l)n,n^l{l)(t){oo))(fx 



3271^ V^^:^ (5 log(/) + log (2/o)) 



3n {ii? — 16) e 



(3.3i; 



b) Contribution of lens-like regions 

Near a; ~ 1 the r.h.c. of eq. (3.30) behaves 3jS 371, I *^ 

1| and the contribution 

of the lens-like regions near x ~ 1, which should be subtracted from the r.h.s. of eq. 
(3.31) is (see (D.9)) 

3.(16 

c) Contributions of the regions Di^, Di^i and D^oo 
We will treat contributions of Di^q and Di^i separately, 
i) Di Q contribution. 

The relevant OPE is 

(p{x)(f)n,n+i{0) = C'{"'3"+^_J^^^)(a;x)-^(0„,„+i + ■ ■ ■ ) 

+ C^{r,3Kni+l)(^^)^"--^-^-"^^-^ (0n,n-l(O) + ■ ■ ■ ) 

It follows from (1.11) that in this region with sufficient accuracy 

Stt / 1 (xx)'^"'"+i~^"'"- 

I{x) ^ — ^ — 

ne \ n + 4: n — A 
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From (A. 5) 

^(n,n+i) ^ _ Vn'^ - 1 {n + 2)(1 - {n - l)e) ^ 

*-'(l,3)(n,n+l)*-^(1.3)(n,n+l)(n,n-l) - ^ H U[e ) 

r(^,n-i) ^ _ Vn^ -l{n-2){l + {n + l)e) , 

<-'(l,3)(n,n+l)<-'(l>3){n,n-l)(n,n-l) - ^^^2 ^ ^'^"^'^^ 

and for the Di q contribution we get 

/(x)(0(x)0(O)0(l)0(oo))c/2x 



-D;,o\-D!„,o 



47r2(n + 2)Vn^^{{n - 8)(1 - (n - l)e) + A{n - 2)elog/) 
3(^2 - 16) (n - 2)n2e2 



^Tl.Ti— 1 + 1 



(3.34) 



_47r2(n - 2)Vn'^ - l((n + 8)(1 + (n + l)e) + 4(n + 2)elog/) 
3(^2 - 16)n2(n + 2)e2 

with two terms corresponding to two intermediate primaries 0n±i- 
ii) contribution. 
The relevant OPE is 

0(a:)0„,„-i(l) = Cl^;;)lZ_,^{{x - l){x - 1))^— + 

+C'{r;KnLl)((^ - - 1))"'' (0n,„-l(l) + 

Considering a; — )■ 1 hmit of (1.11) we get 

?T.e\n-4 n + 4 

Combinations of structure constants relevant for this case are the same as those in 
(3.33). For the A,i contribution we get 

[ /(x)(0(x)0(O)0(l)0(oo))c/2x 

47r2(n + 2)V^2 _ _ 8)(i _ _ i)^) + _ 2)elog/) 
~ 3(^2 - 16) (n - 2)n2e2 



47r2(n-2)v/^^^^((n + 8)(l + (n + l)e) + 4(n + 2)e log(O) 

3(^2 - 16)n2(n + 2)e2 ^ ' 

Again the two terms correspond to two intermediate primaries 0n±i- Notice that due 
to some subtle interplay among quantities involved, for the contribution of Di^i we got 
exactly the same result as for Di q. 
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iii) -Dj^oo contribution. 

Since the structure constant C(i,i)(„,ri-i)(n,n+i) = only the field (j) appearing in u- 
channel OPE gives a nonzero contribution, which is proportional to 

M) r - - 1 (1 - e) , ^,2. ^(.^ 

'-^{l,3)(l,3)*-^{1.3)(n,n-l)(n,ri+l) " 7^ ^ U[e ) (^d.dbj 

Approximating I{x) by I{x) ~ ^ jxl"^"^ we get 

/(x)(0(x)0„,„+i(O)0„,„+i(l)0(oo))c/^x 



,00 

647rVn2 - 1(1 - e(l - 2 log /)) 
3n {p? — 16) 

It remains to collect all the contributions together to get 

167rVri2 - 1 (5n2 - 44 + (n^ + 20) e; 
3n (^2 - 16) (n2 - 4) e"^ 

3.2.3 The matrix of anomalous dimensions 



(3.37) 



+ 0(e°) (3.38) 



There is no need to calculate the remaining two point functions (0n,n+i(l)0n,n-i(O))A 
and (0„,n-i(l)0n,n-i(O))A siucc the former is identical with (0n,n-i(l)0n,n+i(O))A and 
the latter can be obtained from (0n^„+i(l)0„^„+i(O))A by simply replacing n — —n. For 
simplicity of notation let us denote 0n,ri+i = 0i and 0„,n-i = 02, then the two-point 
functions can be represented as 

G'a,/3(x,A) = (0,(a;)0^(O))A (3.39) 



)-^"-^^ - ACS 



The first order coefficients are given by 



= y (0.(i)0/3(o)0(x))c?2x 

r 7r7(6 + A^- A;,)7(6 + Aff- 

= <-{l,3)(a)(/3) ^^^^^ l3.4Uj 

From eq. (A. 2) for the dimensions we have 

A2 = A„,„_, = i + - ^) ^ + ^ K - 1) + O(e^) (3.41) 
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Explicitly, up to 0(e) terms we get 

a, ^ x(n + 2)(2-(2n-l).) ^ ^a, _ x(n - 2)(2 + (2n + 1).) ^ 

Cg = C« = -^^fi3f±^ + 0(e) (3.42) 

and for the second order coefficients we have (see (3.28), (3.38) ) 

(2) _ TT^ (3^3 + 24^2 + 64n + 44) Att'^ {n + 1) {n^ + 7 n'^ + Un + 5) q 



3n(n + 2)2e2 3n{n + 2)'^e 

(2) _ (3^3 - 24^2 + 64n - 44) 47r2(n - 1) (n^ - Tn^ + 14n - 5) q 
^^•2 " 3n{n - 2)2e2 + 3n(n - 2)^6 + > 



ri^)_M2)_ IGyr^y/;?^ (5^2 - 44 + (n2 + 20) 6) 

^^'2 - ^^'^ - 3n (n2 - 16) {n^ - 4) 62 + ^ ^^'^^^ 

Obviously the correlation function (3.39) satisfies the Callan-Symanzik equation 

{xd, + A„ + - eA^A) (x. A) = (3.44) 

As in Section 3.1 let us introduce renormalized fields 

0(f) = i?«,;3(A)0^ 

requiring that the two point functions G^^^ ^^''(0))a satisfy the normal- 

ization condition 

= Se.,^ (3.45) 

In matrix notations we may write 

G^^\x) = B ■ G{x) ■ (3.46) 

Comparing with (3.44) we see that the renormalized two-point function satisfies the 
equation 

2 

(xd, - P{g)d,)G'fl + 5^(r„,,G^) + r^,,GS) = (3.47) 

p=i 

where the /3 function and the renormalized coupling g have been introduced in 
Section 2 and the matrix of anomalous dimensions F is defined as 

r = BAB-^ - eXBdxB-^ (3.48) 
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where 



A 



Ai 
A2 



(3.49) 



Expanding the matrix B up to second order in A 

B = l + XBi + \^B2 + 0{\^) 

5-^ = 1- \Bi + \^{Bl - B2) + 0(A^ 



(3.50) 



imposing the normahzation condition (3.45) and requiring that the matrix of the 
anomalous dimensions (3.48) be symmetric, we find 



2 2e L 



B. 



4 8e 



+ 



A,CW 



1 



A,C(2) 
A,(CW)' 



+ 



A,CW 



1 2 



(3.51) 



Now, all the ingredients to calculate the matrix of anomalous dimensions (3.48) are at 
our disposal. Taking also into account the X-g relation (2.18), we get 

7r^(n + 2)(2-(2n-l)e) vrV 



ri,i = Ai - 
ri.2 = r2,i 



4^3^ 
7r^Vn2 - 1 (2 + e) 



2V3 



n 



7r(7(n-2)(2 + (2n + l)e) Ti^g^ 

i 2,2 — ^2 H \ — 7^ r 



Notice that all matrix elements are regular at e = 0, all double and single poles in e 
disappeared. At the fixed point g = g* (see (2.21)) 



1,1 

r{9*) 
1,2 



2,2 



1 
4 



(2n^ — n — 4) e ^ (n^ — 4n^ + n + 



I6n 



(9*) 



2,1 



- le(l + e) 
2r2 



1 {2n^ + n - 4) e (ra^ + 4n2 + n - 
4 8n 16n 



(3.53) 



The eigenvalues of this matrix 



(3,54) 
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up to O(e^) terms coincide with the dimensions A^^^|j and /S}^_l'j^ of the IR CFT Mp_i. 
We can easily identify also the corresponding normalized eigenvectors and establish an 
explicit map 



(p-i) 



n+l,n 

/,(P-1) 
n— l,n 



1 



(9*) ^ 

n ^ n 



'if 



n 



1 



1 



n n 



(3.55) 



Remarkably the coefficients in (3.55) did not receive neither e nor corrections. Thus 
it is quite perceivable that under the renormalization scheme (3.45), which we have 
adopted following A. Zamolodchikov, the relation (3.55) is exact. The same phe- 
nomenon we will encounter in the next section where we consider a more involved 
case of mixing of the three fields 0n,n±2 and dd(f)n,n- 

3.3 Renormalization of the fields 4'n,n+2i dd'Pn,n and 4'n,n-2 

The OPE 0i,30n,n+2 includes fields from the conformal families [0,i,„+4], [0n,n+2] and 
[(pn,n]- Similarly the product 0i,30n,n-2 produces fields from the families [0n,n-4]5 [0n,n-2] 
and [(pn,n]- Since the dimensions of the primary fields 0n,n±2 and the descendant field 
dd(j)n,n are close to 1 in large p limit, we have a situation when these three fields 
effectively get mixed along RG flow^ [1]. To find the matrix of anomalous dimensions 
we have to calculate all two point correlators of these fields. 

3.3.1 Correlation function (0n,n+2(l)0n,n+2(O))A 

a) Contribution of the region Qi^i^^ 
is given by the integral 



I{x){(j){x)(j)n,n+2{0)(f)n,n+2{l)4>{oo))(fx 



(3.56) 



At large p from AGT relation we got (see Appendix C) 



n,n+2 

(O)0n,„+2(1)0(OO)) 



3x^ - 6x3 + 9^,2 _ ^ ^ 



3(2; — l)^a;^ 



8(3 + n) 
^3(1 +n) 



(2x - 1) (2x2 - 2x + 1) 



4(x - 1)2x2 



(3.57) 



(3 + n)(4 + n) . 
18n(n + 1) ^ ' ^ ' 



^The fields 0n,n±4 have larger dimensions ~ 4 and do not get mixed with these three fields. 
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I{x) can be replaced by ^. Performing integration (3.56) we get 

/(x)(0(x)0„,„+2(O)0„,„+2(l)0(oo))rf2x (3.58) 

167r2 (2^2 + 5n + 1) log(/) 167r2(n + l)log(/o) 
3n{n + l)e 3ne 
7r2(66 + 33n + 128(n + l)log(2)) 2n^{2n + l) 7T^{n + 2) 
2^e ^ 3nPe ^ 6n/^e 

b) Contribution of lens-like regions 
Expanding eq. (3.57) near x ~ 1 we get 

2 + n ^2{n^ + n-2){x + x-2) 



3ra|a;-l|4 3n{n + l)\x - 

^ ^n^ + 5ri + 6){ix-iy + {x-lf) ^ 8(n + 1 
3n(n + l)|x — 1|^ 3ri|x — Ip 

So the contribution of the lens-hke regions near x ~ 1 is (see Appendix D) 
7C f 2 + n f 7T 7T ^^_|_ ^'-^^ + n — 2)tt 



e V 3n V2^o 8/ 3n(n + 1) 

27r(n2 + 5n + 6) 8(n + 1) , / \ , , 

+ ^r^^TV^ + 27r log — 3.60 

3n(n + 1) 3?7, 21q J 

c) Contributions of the regions Di q, Di i and Di^^ 

Contributions of Di^ and Di^i are identical so we compute only Di^ part. Here is the 
relevant part of the OPE 

0(a;)0n,n+2(O) = (XX) ^<^("^3"|„^^i+2)'^«>"+2 + A„,„+2 A 

+ — — -^-j (i + — m::;^ — -^-j 

+ ■■■ (3.61) 

where L^, are the left and right Virasoro generators. To proceed with calculations 
let us notice that the effect of the Virasoro generator L_i or L_i in the three point 
function is rather simple. Namely the relation 

(L_i0„,„(O)0n,n+2(l)0(oo)) = (A„,„ + A„,„+2 " A) (0„,„(O)0„,„+2(l)0(oo)) (3.62) 

and a similar relation with L_i replaced by L_i hold. We see from (1.14) that in this 
region /(x) can be approximated by 

/(x) ~ — |-vr(x + x — log |x|^) — 27re(x + x) log |x|^ (3.63) 
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We need also the combination of structure constants 

4(n + 3)2(1 - (n + 2)e) 



'-"{l,3)(n,n+2)(n,n+2) 



3(n + l)2 



f~i{n,n) ^ _ n + 2 

'-'(l,3)(n,n+2)'-'(".«)("."+2)(l,3) — 

Taking into account that 

A A A r) + 3 

(A + An,n - A„,„,+2) (A„,„ + An,n+2 " A) _ 2(n - 1) (n - + 3)e 



(3.64) 



2A, 



n + 1 



2(n + l) 



we see that the piece of integrand corresponding to the contribution of the family [0„ , 
is equal to 

2 



^ + ^ |^|-4+(3+n)e 

3n 



^ , 2(n-l) (n-l)in + 3)e 



n + 1 



2{n + l) 



(71 
— h TT (x + x — log |x|2) — 27re(x + x) log 



(3.65) 



The part corresponding to the intermediate field 0n,n+2 is simpler. For this one we can 
restrict us with a less accurate expression for I{x) 



I{x) ~ Tflog |x| 



(3.66) 



and the integrand is simply 

4(n + 3)2(1 - (n + 2)e) fn 
3(n + 1)2 

Performing integrations we get 



^ vr log |x|2j |x| 



-2+2e 



/(x)(0(x)0„,„+2(O)0„,„+2(l)0(oo))c?2x 
A,o\Ao,0 

7r2(n + 2)(n + 5)(r2- 1)2 7r2(n + 2) 



3n(n + l)2(n + 3)2f 



3Z2ne 



47r2(n + 2)(n - 1) {{2n^ + 10^2 + 6n - 18) log(/) - - 9n^ - 23n + 1) 

3n(n + l)2(n + 3)2e 
87r2(n + 3)2 87r2(n + 3)2(n + 2 - log/) 



3(n + l)2e2 



3(n + l)2e 



(3.67) 



where the second and the third lines come from the family [0n,n] and the last line, from 
the ^n,n+2 field of the OPE (3.61). 
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Also in this case the contribution coming from the region Di ao is similar to the case 
discussed in Section 2. We should simply take into account that the contribution of 
the field appearing in u-channel OPE is proportional to 



'-^(l,3)(l,3)'-^(1.3)(n,n+2)(n,n+2) 

The result (c.f. eq. (2.10)) is 



4(n + 3)(2- (n + 5)e) 
3(n + l) 



+ 0{e' 



l{x){(j){x)(f)n,n+l{0)(j)n,n+l{^)(p{oo))(fx 



71^ 47r2(n + 3) /2 n + 5-41og/ 
~7P ^ 3(n + l) ~e 

Finally, (3.58) minus (3.60) plus twice (3.67) and plus (3.69) gives 

Stt^ (3n^ + 33n3 + 121n2 + U3n + 20) 

3n(n + l)(n + 3)2e2 
_47r2(n + 5) (5n^ + A5n^ + 143n2 + ISlri + 8) q 
3n(n + l)(n + 3)2e ^ ^ 

3.3.2 Correlation function (0n,n(l)0n,n+2(O))A 

a) Contribution of the region Qij,^ 
is given by the integral 

I{x){(j){x)(f)n,n+2{0)(j)n,n{l)4>{oo))(fx 



The large p limit of the four-point function is very simple (see Appendix C) 

n,n+2 

(O)0„,„(l)0(oo)) 

I{x) can be replaced by 



4 /n + 2 



|x|-2 + 0(e) 



n 



I{x) 



2Tiin + l)e 



+ 5 

and the result of integration is 



4x 



(n + 1)(1 - x) 



2n(n - 3)e 
+ 1 



X 



1 — X 



I{x){(f){x)(j)n,n+2{0)(f)n,n{^)(p{oo))(fx 



l,lQ 



levr^e n + 2 



3(?7, + 5) V n 



((3n-5) log(/) + (n + l)log(2/o)) 



(3.68) 



(3.69) 



(3.70) 



(3.71) 



(3.72) 



(3.73) 



(3.74) 
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b) Contribution of lens-like regions 

We see from (3.72), (3.73) near a; ~ 1 the integrand in eq. (3.71) behaves as 



lne(n + 1] 



n+2 



3(n + 5)|x- 1|2 



So, according to the Appendix D, the contribution of the lens-hke regions near a; ~ 1, 
which should be subtracted from the r.h.s. of eq. (3.74), is 



16.^(n + l). In^2^fi ^^^^^ 



3(n + 5) V n \2lo 

c) Contributions of the regions Di^, Di^i and D^oo 

Let us compute contributions of D^o and Di^i separately. 

i) Di Q contribution. 

The relevant OPE already appeared in (3.61). Instead of (3.62) now we need the 
analogous relation 

(L_i0„,„(O)0n,„(l)0(cX))) = {2An,n - A) (0n,n(O)0n,n(l)0(oo)) (3.76) 

For I{x) the first equality in (1.11) should be used. When calculating the contribution 
of (j)n,n+2 the hypergeometric functions can be simply replaced by 1, but for the con- 
tribution of the family [(j)n,n] they should be expanded up to linear orders in x or x. 
Below we present the expressions of relevant combinations of the structure constants 
with required accuracy 



*-^(l,3)(n,ri+2)*-^(l,3)(n,n)(n,n) - y — ^ H U[e ) 

Mn,n+2) ^ ./^ (^ + 3)(2-(n + 2)6) , 
<-(l,3)(n,n+2)'-(l,3)(n,n+2)(n,n) = ^ — S^^T+Y) ^ ' ^ ' 

For the Di o contribution we finally get 

I{x){(j){x)(f)n,n+2{0)(pn,n{^)<Pioo))(fx 



L 



Di,o\DiQ,o 

^2/^ _ p . /n+2 



n 



3(n + 3)(n + 5) \ \ 2 



1 + e ( n + - + (n + 3) log / 



3{n + 5) V V2 



n 



1 + e -+4 + 21og/ (3.78) 
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with two lines corresponding to the [0n,n] and (pn,n+2 contributions respectively. 

ii) Di i contribution. 
The relevant OPE is 

</'(x)0„,„(l) = \X- ir(^---^-'-^)C^{r5S)</'n."+2^ (3.79) 
+C^(,3)(.„)I- - 1| (1 + j (1 + J 0n.(l) + ■ ■ ■ 

The relation between the three point functions relevant for this case is 

(L_i0„_„(l)0„_„+2(O)0(oo)) = (A - A„,„ - A„,„+2) (0„,„(l)0„,„+2(l)0(oo)) (3.80) 

Notice the flip of sign compared to (3.62) due to interchanging of the points and 
1. Now for I{x) the second equality in (1.11) should be used. The hypergeometric 
functions should be expanded around x = 1. When calculating the contribution of 
(f>n,n+2 it is enough to keep only the constant term but for the contribution of [4>n,n] 
also the terms linear in a; — 1 or x — 1 should be taken into account. During calculation 
one encounters the same combinations of the structure constants as in (3.77). Finally 
we get for the Di^i contribution a result identical to that of Di^ given by (3.34). The 
same phenomenon we encountered earlier in Section 3.2.2. 

iii) -D;,oo contribution 

Only the field appearing in u-channel OPE gives a nonzero contribution. Since 



r _ n + 2 2(2 -3e) , 



and, from the third equality in (1.11) 



JU) ^ _!![!(!^_J^ (l + (n + 5)e) \x\'^' (3.82) 
n + 5 



we get 



/(x)(0(x)0„,„+2(O)0„,n(l)0(oo))rf^a; 

-Di,oo\£'!Q,oo 

n + 2 87r2(n-3) (2 + (2n + 7 + 41og/) e) 
n 3(n + 5) 

It remains to combine all the contributions. The result is 



(3.83) 



n + 2 47i\n- 1) (6n + 22 + (5rz^ + 37n + 64) e) , 
n 3(n + 3)(n + 5) 
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3.3.3 Correlation function (0n,n+2(l)0n,n-2(O))A 
a) Contribution of the region Qi^i^^ 



l,ln 



/(a;)(0(a;)0„,„_2(O)0„,„+2(l)0(oo))d^x (3.85) 
The large p limit of the four-point function is (see Appendix C) 



{(f>{x)<f>n,n-2mn,n+2im00)) = ^^^^ + 0(e) (3.86) 



I{x) can be replaced by 



and the result of integration is 

/(x)(0(x)0„,„_2(O)0n,„+2(l)0(oo))(i^X 



^64 log (2/ /o) + 33-1-4) (3-88) 



b) Contribution of lens-like regions 

From (3.86) and (3.87), near a; ~ 1, the integrand in eq. (3.85) behaves as 

4n (1 - 2(x + X - 2) + (x - 1)^ + (x - 1)^ + 4|x - Ip) 
?)ne\/v? — 4 |x — 

Consequently, from the Appendix D, we see that the contribution of the lens-like regions 
near x ~ 1 is 

4vr f TT TT TT , s. , s. , f I W 

2(vr) + (27r) + 87rlog — (3.89) 



c) Contributions of the regions Di q, Di i and Di ^o 
i) Di Q contribution 
The relevant OPE: 

0(x)0„,„-2(o) = (xx)^--^---^dr;3l,„_2) X (3-90) 



' + 2A 



^-') + ^ ^ ~ "^^"0 ^ 



-27- 



Impact of L_i on the three-point function: 

(L_i0„,„(O)0„,„+2(l)0(oo)) = (A„_„ + A„,„+2 - A) (0„,„(O)0„,„+2(1)0(OO)) (3.91) 

In the first expression of (1.11) for I{x), the hypergeometric functions should be ex- 
panded up to linear orders in x or x. 
Relevant combination of the structure constants: 



^(l,3)(n,n-2)'-^(l,3)(n,n)(n,n+2) - ^ H ^[(^ ) 



Final result for the Di^q contribution: 

I{x){(f){x)(l)n,n-2iO)(pn,n+2{^)<P{oo))(fx 



(3.92) 



A,o\Ao,o 



W flO n2 + l + 2(n2-9)log/\ 

(3.93) 



3n (^2 - 9) -4 V 4/2e 



ii) contribution 
The relevant OPE: 

0(x)</.„,„+2(l) = k - ir(^"--^-"+^-^)C^(t3K.,n+2) (3-94) 

X f 1 + (X - 1)L_, 

X ( 1 + 1^ ■ - l)L_i J (PnA^) + ■■■ 

Impact of L_i on the three-point function: 

(L_i0„,„(l)0„_„_2(O)0(cx))) = (A - A„,„ - A„,„_2) (0„,„(l)0„,„_2(l)0(oo)) (3.95) 

The combination of structure constants required for this computation coincides with 
that given by eq. (3.92). The explicit calculation shows that in this case too, the Di i 
contribution is identical to that of Di^q, given by (3.93). 

Note also that the contribution of Di^^o is negligible. Combining all the contributions 
for the case at hand we get 

Oi^) (3.96) 



3e2n(n2 - 9)V^2 _4 
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3.3.4 The matrix of anomalous dimensions 

The remaining two point functions (0n,n-2(l)0n,n-2(O))A and (0„,n(l)0n,n-2(O))A can 
be obtained from (0„,n+2(l)0n,n+2(O))A and (0„,n(l)0n,n+2(O))A replacing n by -n. So, 
we have all necessary material to repeat the steps of Section 3.2.3 and calculate the 
matrix of anomalous dimensions for the fields 

4>1 = 4>n,n+2 

02 = (2A„,„(2A„,„ + 1))-' dd(j)n,n 

03 = 0n,n-2 

The two-point functions of these fields can be represented as in (3.39), but the indices 
now take the values a, (3 = 1,2,3. Replacement of the field 0„_„ by 02 in a two point 
function, at a given order k of the perturbation theory, results in an extra multiplier 
which is easy to calculate. Here is the rule: the coefficients Cj^l = C'^'l-, for a 7^ 2, and 
C2 2 should be endowed with extra multipliers 

(fee - A, - A2)' 
2A„,„(2A„,„ + 1) 

and 

/ (fc6-2A2) (A:6-2A2-l) y 
V 2A„,„(2A„,„ + 1) ) 

respectively. The numerators come from the derivatives and the denominators from 
the normalization factor, present in the definition of the field 02. The dimensions at 
the zero coupling A = are 

Ai = A„,„+2 = 1 - e + — — + 0{^) 

Z lb 

- I 

A2 = 1 + A„,„ = 1 + — — + 0(e=^) 

lo 

A3 = A„,„_2 = 1 + e + + 0(6=^) (3.97) 

/ lb 
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For the first order coefficients we get 



3n 

{n+T){n + 3)e 



V3(n^ — l)e 



Ml) _ Ml) _ V 3n s 

Ci-)^^-(--y^ + (-^^-^0(e) (3.98) 
y3(n - l)e ^ ' ^ ' 

From (3.70), (3.84), (3.10), (3.96) and the above presented considerations, for the 
second order coefficients C^^^ = C'p'\ we find 



rt-2 



(2-6) 



^(2) 



'-"1,2 



^(2) 



c. 



(2) 
2,2 



^(2) 
'-'2,3 



c. 



(2) 
3,3 



877^ (3ra^ + 33^3 + 121^2 + 143n + 20) 

3n(n + l)(n + 3)2e2 
47r2(n + 5) (5n^ + Abn^ + 143n2 + 151n + 8) 



3n(n + l)(n + 3)2e 



+ 0(e°) 



647r^ 



n+2 



+ 



(3ra + ll) 327r2j^(n2 + 18n + 57) 



3(n + l)(n -f 3)(n + 5) e2 3(n + l){n + 3)(n + 5)e 

3207r2 3207r2 „ 

; , + 0(e") 

3n(n2-9)Vn2-4e2 3n(n2 - 9) _ 4 ^ 



+ 0(e°) 



1287r2 167r2 (n^ + 19) 

3 (n2 - 1) e2 3(^2 - 1) e 

647r%/^(3n-ll) 



-f 0(e°) 



327r\/^ (n2-18n + 57) 



+ 



3(n - l)(n - 3)(n - 5)e2 3(n - l){n - 3)(n - 5)( 
87r2 (3n^ - 33^3 + 121^2 - 143n + 20) 

3n(n- l)(n-3)2e2 
47r2(n - 5) (5n^ - 45^^ + 143^2 - 151n + 8) 



+ 0(e°) 



3n(n — l)(n — 3)2e 



+ 0(6°) 



(3.99) 
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With this input we can repeat the procedure we carried out in Section 3.2.3 and compute 
the matrix of anomalous dimensions. Here is the final result: 

^ ^ , ng{n + 3)i2-in + 2)e) , 87rV(n + 2) 



ri,2 


= r2,i ^ 


ri,3 


= r3,i ^ 


^2,2 


^ A2 + - 



^9{n - 1)^ (2 - e) An^g^n - l)y^ 



n+1 ^ 3(n + l) 

27r^(4 - {n^ + l)e) 47r2^2 (^2 ^ 3) 
V^(n2- 1) ^ 3(^2 - 1) 

7ro(n-3)(2 + (n-2)e) 87rV(^-2) 

r3,3 ^ A3 + ^ \- ^ + J ^ 3.100 

V3(n-1) 3(n-l) 

Again we see that all matrix elements are regular at e = 0, all double and single poles 
in e disappeared. At the fixed point g = g* (see (2.21)) 

~ ^ 2{n + l) + 16(n + l) ^ 

, , (n-l).f^ e(e + l) 



Here are the eigenvalues of this matrix 



a(^*) 


= 1 + 


{n + 1} 
2 


^2 


= 1 + 


{n^ - 1 
16 


^3 


= 1 - 


(n-l) 
2 



(„+!)(„ + 7).' 

16 ^ ^ 

^2 

.3n 



- + 0(e^ 

+ ^^^^^4p^ + 0(e3) (3.102) 
16 

which, up to 0{e'^) terms coincide with the dimensions A^^2^!i5 1 + ^n,n^'^ and A^^_2\ 
of the IR CFT Mp_i. It is easy to find the orthogonal matrix which diagonalizes the 
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matrix of anomalous dimensions (3.101) and to establish the explicit map 



2 ln+2 

' + n n + 1 



n + l n^ — 1 72 — 1 



'■'n-2,n — Vl T- V2 ' , 

n n — 1 [n — Ijn 



+ 7 TV- n (3.103) 



where we have introduced the notation 

cf>t'^ ^ (2A(^;i)(2A(^;i) + 1))-' 990^^) (3.104) 
In this case too we see that the coefficients in (3.103) do not receive e or corrections. 
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A Minimal models 

For the readers convenience we present here few facts about unitary series of the min- 
imal models [5, 8] denoted by Mp, p = 3,4, . . . . The central is given by 



6 

p{p+l) 



This theory has finitely many spin-less^ primary fields denoted by 0„,m, where n G 
{1,2,..., p—1}, m E {1, 2, . . . , p} with conformal dimensions (the famous Kac spectrum 
[6]) 

A_ = ^^^^ + ^-^?^ (A.2) 
4 4p 4(p+l) ^ ^ 



^We consider here the so called diagonal series only. 
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There is an identification (f)p-n,p+i-m = 4>n,m so that the number of primary fields is 
equal to p{p — l)/2. The field with dimension is the identity operator. The 
operator product expansions satisfy the fusion rules 

ni+n2 — l mi+m2 — 1 

0ni,mi0n2,m3 ^ [^na.ms] (^-3) 

"■3 = |'^i— "-21+1 m3=|mi— m2|+l 

ni+n2—n3S:2'Z+l mi+rri2— m.3£2Z+l 

In this paper we consider perturbation of the minimal model Mp by the relevant field 
0i_3. It's dimension 

Ai,3 = 1 - = 1 - e < 1 (A.4) 

p+l 

For large p this field becomes nearly marginal which is the main reason why in this 
region the non-trivial RG behavior can be investigated by means of the perturbation 
theory. The structure constants of the OPE have been computed in [9]. A slightly 
more compact expression which we present below is taken from [10] 



X 



'7 (p - l)7(mi - nip ^)7(m2 - ri2P"^)7(-m3 + ngp-^) 
7(1 - p-i)7(-ni + mip)7(-n2 + m2p)'^{n^ - m^p) 
t 

n n ~ ^P^^'^ + ^^3 - (i + 'm^)p)ii -ni- (j - mi)p){i - ^2 - (j - m2)p))" 



s t 

X 

i=l j=l 



xjj7(zp ^)7(-m3 + (i + n3)p ^)7(mi + (i - ni)p ^)'^{m2 + {i - n2)p ^) 
1=1 
t 

X YllUphi-ns + U + ms)p)'y{ni + (j - mi)p)7(n2 + (j - "^2)^) 
j=i 
where 

/ , r(a;) 77,1+^2-^3-1 , mi + m2 - m3 - 1 

^^^^"f(r^' ' = 2 ' ' = 2 

B Computation of I{x) 

One way to get the result (1.10) for the integral (1.8) is to notice that I{x) satisfies the 
hypergeometric differential equation independently with respect to the both variables 
X and X. The starting point is the identity 

(a;(l - x)dl + {l-a-c+{a + b + 2c- 2)x)d^ + c(l — a — 6 — c)) 

X (y"-i(l - yf-'iy - xY) = dy {cy\l - y)\y - xf-') (B.l) 
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which shows that as a function of the variable x, I{x) is a hnear combination of the 
hypergeometric functions 

F(l — a — b — c, — c, 1 — a — c, x) 

and 

x''^''F{a,l-b,l + a + c,x) 

The same conclusion is true also for the conjugate variable x. The condition that 
the function I{x) is single valued around the points x = a x = 1 fixes the correct 
combination of holomorphic and anti-holomorphic parts up to a constant which in its 
turn can be easily evaluated considering the special case a; = 0. The final result is 
presented in eq. (1.10). 

C Four-point functions at large p limit 

Since the structure constants of OPE for the minimal models are known (see A. 5), 
to construct the correlation functions it remains to calculate related conformal blocks. 
According to AGT relation [11] this conformal blocks in a simple fashion are related 
to the instanton part of the Nekrasov partition function of = 2 SYM theory with 
SU (2) gauge group with four fundamental hypermultiplets. In the large p limit the 
minimal models approach to a free theory (the central charge c ~ 1), so it is not a big 
surprise that in this limit conformal blocs of degenerated primary fields are very simple 
and can be expressed in terms of rational (and also logarithmic in the cases when 
leading corrections in 1/p is required to be taken into account) functions of the the 
cross ratio of the coordinates. It is straightforward to compute Nekrasov partition [12] 
function up to desired order in instanton expansion using combinatorial formula found 
in [13] and extended to the case with extra hypermultiplets in [14]. Computing first 
few coefficients of the instanton expansion (for more confidence we made calculations 
up to 6th order ), adjusting appropriately the parameters in order to get the required 
conformal block and taking the large p limit one can easily guess the exact dependence 
of the conformal block on the cross ratio of insertion points (which is the same as the 
instanton counting parameter, from the gauge theory point of view). In this way we 
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got expressions for the correlation functions 



(01,3(a;)01,3(O)0n,n.(l)0n,n.(oo)) = \x 



4e-4 



{n' - 1) 
12|x|4 



X 



+ 



2 1-x 2 1-x 



log'll-xM +0{€') 



(01,3(a;)01,3(O)0„,n+l(l)0n,n+l(oo)) 



1 — X 



(01,3(x)0i,3(O)0„,„_i(l)0„,„+i(oo)) 

(01,3(x)0i,3(O)0„,„+2(l)0n,n+2(oo)) 

2 



— x) 
4Vn2 - 1 



+ 
1 - 



2(n + 2) 



3?2 
2 



1 - - 
2 



x(l — x) 



+ 0{e) 



xil — x) 



+ 0{e) 



1 - 2x + - 2a;3 ^ ^ 



x2(l — x)2 



+ 



8(n + 3) 
3(n + 1) 



3x 



_|_ i£l 

2 -T J- 4 



x(l — x)2 



(n + 3)(n + 4) 
18n(n + 1) 



X 



1 — X 



+ 0(6) 



(01,3(x)0i,3(O)0„,n(l)0n,n+2(oo)) = 
(01,3(x)0i,3(O)0„,n+2(l)0n,n-2(oo)) 



n 



3 V ^ 
\/rfi - 4 
3n 



X 



X 



1 — X 



0(e) 
' + 0(e) 



(C.l) 



As a nontrivial check, we have tested the crossing invariance of all these correlation 
functions. Of course the interested reader can get convinced in correctness of our 
expressions also by examining the third order differential equation satisfied by any 
conformal bloc which includes the degenerated field 0i 3 [5]. 

Performing the conformal map x — )■ 1/x with the help of the eq. (1-6) we get the 
correlation functions (2.1), (3.2), (3.20), (3.2), (3.30), (3.57), (3.72), (3.86) used in the 
main text. 



D Integrations over lens-like regions 

Here we compute the contributions of lens-like regions 

Dl = 1^1-/0,0 n A,i 
Dr = /^i+/o,o n Di^i 

(D.l) 

(see Fig.l and also the discussion coming after the eq. (2.4)). Using Green's theorem 
the integrals over lens-like regions can be easily transformed to the contour integrals 
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over their boundaries. The integrals over the arcs which belong to the circle |a; — 1| = / 
are trivial. Instead, the integrals along remaining parts of the boundary which lay on 
|x| = 1 + /o (for the right lens-like region D/j) or on = 1 — Iq (for the left lens-like 
region Dl) seem more complicated, but fortunately these contour integrals too (with 
an exception to be considered later) admit exact treatment. Below we give the details 
on the integration along the arc |a;| = I + Iq. The formulae for the other arc |x| = 1 — /q 
can be found by a simple replacement Iq ^ —Iq. During calculations we heavily employ 
the formulae (we use the notation r = \x — l\ and the angles 0, a are depicted in Fig. 2) 




(1 + /o) sin(a) = sin(y9); r = (1 + Zq) cos(a) — cos(v9) (D.2) 
(1 + Iq) cos(a) + cos{(p) 



loik + 2) 

First note that for the arbitrary region I) 

(Px f drdcp 1 f dip 



V |a; — 1|^ jgj) I jgj) 



V 



(X-1)2 (X-I)V Jgjy r 

d^x (I I \ f 2dsm{(p) 



(D.3) 



1| \x 1 X 1/ J 

(D.4) 

When restricted on the circle = l + /o the one forms on the r.h.s. can be represented 
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as total derivatives. Thus 



d^p (1 + IqY cos^(q;) + cos^((/?) + 2(1 + Zq) cos(q!) cos((^) 

/^(l + /o)2 

_ (1 + /o)^d(a + + (1 + /o)c?sin(a + Lp) 
sin(29?)rfr (1 + kf 



d(p 



r 2 



d{2a - sin(2Q;)) 



r V ^0^0 + 2) ; ^ ' ' 

With above formulae at hand it is easy to evaluate the integrals over lens-hke regions 
in the limit Iq/I — > and I — > 

+ ^ (D-6) 



d^x ( \ + \ 1 27r 

r d^x ( \ 1 

JdlUDr 



+ 7 ~ TT 



X — IrVa; — 1 X — 1 



(D.7) 



We need also the integral 



^Xi P ~ -"-I ^sx> 



Unlike the previous cases this integral can not be evaluated exactly in terms of elemen- 
tary functions. Nevertheless it is not difficult to show that up to terms vanishing in 
the limit Iq/I ^ and Z — > it is equal to 



27rlogf4) (D-9) 
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